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ABSTRACT 

Consideration is given to the fully developed heat transfer 
for longitudinal laminar flow between cylinders arranged in an equilateral tri- 
angular array. The analysis is carried out for the condition of uniform heat 
transfer per unit length. Solutions are obtained for the temperature distribution, 
and from these, Nusselt numbers are derived for a wide range of spacing-to-diameter 
ratios. It is found that as the spacing ratio increases, so also does the wall-to- 
bulk temperature difference for a fixed heat transfer per unit length. Correspond- 
ing to a uniform surface temperature around the circumference of a cylinder, the 
circumferential variation of the local, heat flux is computed. For spacing ratios 
of 1. 5 ~ 2. 0 and greater, uniform peripheral wall temperature and uniform periph- 
eral heat flux are simultaneously achieved. A simplified analysis which neglects 
circumferential variations is also carried out, and the results are compared with 
those from the more exact formulation. » 

NOMENCLATURE 
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A,B, ...,F constants in the temperature solution 

Cp specific heat at constant pressure 
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d e equivalent diameter, eq. (27) 
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mass flow parameter, (w/l2py s^ 
summation index 
normal direction 
Nusselt number, hd/k and hd e /k respectively 
summation index 
static pressure 

heat transfer rate per unit length from one cylinder 
local heat transfer rate per unit area 
circumferential average of q, q/2jtt 0 

radial coordinate; r Q , cylinder radius ; r*, outer radius of 
equivalent annulus 

half spacing between centers of cylinders 
temperature; T v , wall temperature, T^, bulk temperature; 

T 0 , entering fluid temperature 
axial velocity distribution 

mass rate of flow associated with one cylinder 
axial coordinate measuring distance from the entrance of the 
passage 

dimensionless coefficients in velocity solution 
angular coordinate 

i 

absolute viscosity 
density 

dimensionless coefficients in temperature solution 


1 dp\ 

H dz/ 


INTRODUCTION 

In a previous investigation (ref. l), consideration was given to the Ion 
gitudinal laminar flow between cylinders arranged in regular array. Here, at 
tention is directed to analyzing the heat transfer characteristics of such a 
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flow. This configuration has potential application in compact heat exchangers 
and has heen frequently considered in connection with the cores of nuclear 
reactors. 

A schematic diagram of the system under study is shown in Fig- 1. The cyl- 
indrical rods between which the fluid flows are arranged in equilateral tri- 
angular array, and the flow is laminar and full-developed. The thermal situation 
to be considered here is the condition of uniform heat transfer along the length 
of the passage. The goal of this analysis is to determine the fully-developed 
heat transfer characteristics of the problem, taking into account the circumferential 
variations of the temperature distribution which inevitably arise in the cross 
section of a noncircular passage. The existance of such peripheral variations 
is what makes the problem both interesting and challenging. The starting point 
of the study is the basic law of energy conservation appropriate to the three- 
dimensional situation. Solutions are found for a wide range of spacing to diam- 
eter ratios, and corresponding heat transfer results are obtained. For large 
spacings, the flow about any given cylinder should be little influenced by the 
presence of neighboring cylinders. With this in mind, a simplified analysis has 
been carried out neglecting circumferential variations, and the heat transfer 
results thus obtained are compared with those from the more complete formulation. 


ANALYSIS 

Mathematical formulation . - To predict the heat transfer characteristics 
from first principles, we must begin with the fundamental physical law of energy 
conservation. The mathematical representation of this law appropriate to a 
fully developed flow may be written in cylindrical coordinates as follows 

(i) 


St 

pc p u ^ 


_ fb 2 T 1 

= W 2 + r 


St _i_ S^t S 2 t 
^ + r 2 be 2 + Sz 2 
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where u, the axial velocity, depends on both cross section coordinates r 


and 0. 



4 


To facilitate applying Eq. (l) to the system under study, we refer to 
Fig. 2. The left hand sketch is a cross-sectional view. From the symmetry of 
the situation, it is easy to see that consideration need he given only to the 
typical 30° element which is blackened in the sketch. Once the solution is 
known within this element, then it is known everywhere throughout the cross 
section. An enlarged view of the typical element is shown at the right, along 
with the boundary conditions and dimensional nomenclature. The condition that 
the normal derivative be zero is an expression of symmetry. On the inner bound- 
ary, r » r Q , the temperature is assigned a value T v which is independent of 
angle. Physically speaking, this means that the surface temperature is uniform 
around the circumference in a given cross section, although it will vary along the 
length. Of course, there still remains a circumferential temperature variation 
within the fluid and this will be included in the analysis. Before leaving the 
boundary conditions, it is interesting to observe that the case of uniform pe- 
ripheral surface temperature is one of two interesting limiting conditions. * 

The other, uniform peripheral heat flux, still remains to be analyzed. 

The first step in attacking Eq. (l) is to introduce the condition of fully- 
developed heat transfer. For the case of uniform heat transfer per unit length, 
this condition is. 


= constant - -9- . MW- 

o z wcp (w/12;cp 


( 2 ) 


where w and Q respectively represent the mass flow and the heat transfer per 
unit length associated with a single cylinder in the arrayj while (w/l2) and 
Q/l2 are the quantities appropriate to the typical 30° element of Fig. 2. 

* 

A thorough discussion of thermal boundary conditions in noncircular ducts- 


is given in Ref. 2 
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Next, it is necessary to know the velocity distribution u. In Ref. 1, a 
very accurate, hut approximate solution for the velocity has been derived in the 
form 
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The coefficients 8^ through 67 have been tabulated to four decimal places in 
the reference for a wide range of spacing ratios, and additional significant fig- 
ures were available when needed in the heat transfer analysis. Another piece of 
information which is needed from the velocity solution is the relation between 
the mass flow w and the pressure drop. For convenience, we define a dimen- 
s ionless grouping M as 


M = 


(w/l2p) 



(4) 


Numerical values of M may be obtained as a function of spacing ratio from 
Fig. 6 of Ref. 1. 

Then, the information contained in Eqs. (2) through (4) may be introduced 
into the energy Eq. (l). After rearrangement, the governing equation for T is 
found to be 


d 2 T + 1 dT 1_ c#T 
Sr 2 r ^ r 2 Se 2 


Q/l2k ) -J3 
Ms 2 ) * 




(5) 


Solution of the governing equation. - In approaching Eq. (5), we look for 
particular and homogeneous solutions. Thus, the solution for T is expressed 
as the sum of these separate solutions as follows 

T = T p + Th (6) 


*In Ref. 1, the symbol Q is used to denote w/l2p. 
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For a particular solution, any function satisfying Eq. (5) will suffice. It is 
easily verified by direct substitution that the following expression is a satis 
factory particular 1 solution 



To find the homogeneous equation, we set the right side of Eq. (5) equal to 
zero and get 

1 S 2 T h 

ar^ r de*' 

This is the well-known Laplace equation, the solution to which may he found in 
various hooks on advanced calculus as 

T h = A + B k r + ^ (C m r m + D m r " m )( E m cos + F m sin m0 ) ( 8 ) 

m 

where m takes on integral values to insure that T is single valued, i.e. , 
that T (0) - T(Q + 2*). The constants A, B, . .., F remain to he determined 
from the boundary conditions. 

So, the general solution of Eq. (5) is given hy Eq. (6) in conjunction 
with Tp and T^ from (7) and (8). But, to complete the solution, there 
still remains the task of satisfying the boundary conditions. Starting with 
the most easily applied conditions, we first Impose the symmetry requirement 
that dr/dO = 0 at 0 = 0° and at 6 = 30° (see fig. 2(h)). After differ- 
entiating Eq. (6) and substituting from (7) and (8), it is found that the con- 
dition at 0=0° yields 

F m = 0 (9) 

while from the condition at Q = 30°, it follows that 

m = 6, 12, 18, ■ . . 


( 10 ) 
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T = T p +A+Blnr + ^ (C^r 6 ^ + D^^Jcos 6j0 


to insure that sin mjt/6 = 0. With this information, it is convenient to re- 
phrase the temperature solution (6) as 

( 11 ) 

j=l “ 

where Tp is still given by Eq. (7) and Ej has been set equal to unity with- 
out loss of generality. 

Next, it is required that at the surface of the cylinder, r = r Q , the tem- 
perature take on the value T w independent of angle. Applying this condition 
to Eq. (ll), it is found that 


A = T, 


•w 


B In r. 


(Q/lEk)r| j 3 ( r °Y*! fi-.l 

° ' — s* - r =" + ^ s (12a) 


Dj = - Cj^ - 
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12n+2 


n = j 


(12b) 


(13) 


Ms 2 2(36n 2 - l) s 6n 

Further, to make certain that energy conservation is satisfied on an overall 
"basis for the typical element of Fig. 2(h), we equate the heat transfer per 
unit length, Q/l2, to the heat conducted at the surface of the cylinder, i. e. , 

.jt/ 6 

”0 VS? 'r 0 

This insures that there is no net energy transport across the other boundaries 
of the typical element. Of course, overall energy would be automatically satis- 
fied by an exact solution^ but if an approximate solution is being contem- 
plated, then it must be imposed as a separate requirement. Inserting Eq. (ll) 
into the overall conservation condition (13) gives 


/•V 6 

«/ 12 - / - k (il r ° de 


B 


- F. ^ + & r 2 ^ 41 

Ms 2 L it 4* r o ~ 16 s 2 ^ 


(14) 


Before going ahead to the final boundary condition, it is useful to bring 
together the findings of the previous paragraph. Introducing Eqs. (12) and 
(14) into the temperature solution (ll) and substituting for Tp from Eq. (7) 
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gives 



(15) 


Remembering that 5n and M are known numerical constants whose values depend 
upon s/r 0 , it is seen that the temperature distribution will "be completely de- 
termined as soon as the Cj are found. At our disposal is the condition that 
St/SN =0 on the right hand boundary of the typical element (see fig. 2(b)), 
on which r = s/cos 0. In applying this condition, it is convenient to use the 
identity 


St St n St sin e 

Sif - 5F cos e - Se 


(16) 


Then, introducing the temperature distribution (15) and setting St/Sn = 0 for 
r = s/cos 0, there is obtained after considerable rearrangement 


£«*(- 0) 1_6J [coB(6j - 1)0 + ^ r ° . ° OS . 


e\ 12 J 
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“ - Z | 56 J n . 1 (t) (cos e)6n+1 + 1)8 

cos(6n - l)0 + (l/5n)cos 0 cos 6n0 
4(6n + l) ( cos 0) 6n+1 

(r 0 /s) 12n (cos 0) 6n-1 
4(1 - 6n) 
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cos' 


Q:os(6n + l)0 - (l/3n)cos 6 cos 6n0]j 

e . ifiaf . jH \ + i + a2 

8\ s / 2n 111 [v 0 cos Q) 16 cos 2 0 4jt 


( 17 ) 



where 




p S 6j M 

3 (Q/l2k) 


(17a) 


The only unknowns appearing in Eq. (17) are the coefficients au (i.e . , C<). 

In considering ways of solving for the cjoj ^ the first thought might he to use 
the Fourier series method. If workable, this procedure would provide an in- 
finite number of values while satisfying Eq. (17) at all points along the 

boundary, 0° — Q — 30°. However, Fourier methods cannot be applied in the pres- 
ent problem and some other approach must be used. The procedure employed here 
is as follows: First, the summation involving the o)j is truncated after i 

terms. Next, the condition (17) is applied successively at i values of 6 
in the range 0° — 6 — 30°. This yields i linear equations, each of which 
contains i unknown values of ooj. As soon as the s/r Q ratio is specified, 
the right hand side of each equation can be reduced to a numerical constant 
since the values of ^ and M are available from Ref. 1. So, there are i 
linear, inhomogeneous algebraic equations, the solution to which yields nu- 
merical values for the coj* To check the accuracy of the results, the entire 
procedure may be reapplied using additional terms in the summation. From this, 
it is found that the additional terms do not significantly affect the numerical 
values of the first few coefficients. Another important fact is that the mag- 
nitudes of the ooj decrease quite rapidly with increasing j, so that only the 
first few ooj are really important. These favorable circumstances somewhat 

ease the task of finding a sufficient number of accurate (&. values. The nu- 

J 

merical results which have been obtained for the cd- are listed in table I as 

J 

a function of the spacing ratio s/r 0 . 

With the tDj values now known, then the coefficients Cj are also known 
through Eq. (17a). With these, we may return to the temperature distribution 
(15). For convenience, the summation involving the Cj may be rephrased as 
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I 

j=l 




- ^ i idf t ■ i^r] 

j=i 


cos 6j0 (18) 


Since ^ and M are also known numerical constants (ref. l), then it is 
clear that once the geometrical parameters r Q and s have been specified, 
the temperature distribution within the typical elements is determined. So, we 
are now in a position to proceed with the determination of the heat transfer 
results. 


HEAT TRANSFER CHARACTERISTICS 

Wall-to-bulk temperature difference and Nusselt number . - The solution 
which has just been presented gives the distribution of the temperature in a cross 
section relative to the wall temperature T w in that cross section. So, to 
complete the solution, there still remains the task of finding T w . Furthermore, 
in practice, it is the wall temperature itself which is usually the quantity of 
greatest practical interest, for example, on a materials limit or a thermal 
stress basis. 

Now, in a situation where the wall heat flux is prescribed, the bulk tem- 
perature is immediately determined according to the following relation 

Tb - T 0 + ^3- z (19) 

where T 0 is the temperature of the fluid at the entrance to the heated pas- 
sage. Making use of this, the wall temper ature may he written as 

T w = Tb + (T v - Tb) = T 0 + JL z + (T w - Tb) (20) 

From this, it is clear that the wall-to-bulk temperature difference (T w - T^) 
is the key to the determination of the wall temperature. In addition, as soon 
as T w - Tb is known, then Nusselt number results can he obtained. 

To compute the wall-to-hulk temperature difference, we start with the 
standard definition of the hulk temperature 
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ff 


Tpur dr <10 


T b ■ 


Area 


ff pur dr <3.0 
Area 


( 21 ) 


The denominator is simply the mass flow. Specializing to the typical element 

of Fig. 2(h) and subtracting away the wall temperature, there is obtained 

' s/cos 0 

-o 


- Tb = 


/ 30° „ . . 

/ x (T w - T)ur dr d0 

= 

w/12 


( 22 ) 


In evaluating the integrand, the velocity distribution is introduced from 
Eq. (-3), while the temperature distribution is taken from Eq.. (15) as modified 
by (18). After making these substitutions and rephrasing in terms of dimen- 
sionless variables, the following functional relation is revealed 


- Tb 
Q/12k 


f(B/r c ) 


(25) 


This is an interesting result in that the dependence on geometry enters only 
as the ratio s/r Q , rather than as s and r Q separately. It is easy to see 
that multiplying out the temperature and velocity distributions in the integrand 
of Eq. (22) leads to an exceedingly lengthy expression which Is made even longer 
when the double integration is carried out. To display only the end result of 
the integration would require considerable space. In the interest of a concise 
presentation, the details associated with the evaluation of Eq. (22) will be 
omitted and only a few general remarks will be made as follows: In organizing 
the computations, separate account was kept of the contribution due to the non- 
series portion of u and T - T w and of the contribution due to the series 
terms. As expected, the relative importance of the series terms increased as 
the spacing decreased. The integrations which arose in connection with Eq. (22) 
were carried out analytically as far as possible} however, it was ultimately 
necessary to compute some specific integrals by numerical means and these are 
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given in the Appendix with the hope that this information may he useful to 
others. 

From the evaluation of Eq. (22), numerical results for the wall -to -hulk, 
temperature difference have been obtained as a function of the spacing ratio 
s/r Q . These results will be reported in terms of the customary Nusselt number 
representation. First, a heat transfer coefficient is defined as 


= Q - Q/12 

- rtd(T w - T-b) (rtd/12) (T w - T b ) 


(24) 


Then, using the cylinder diameter as characteristic dimension, the Nusselt num- 
ber becomes 



It may be seen that the Nusselt number is essentially the reciprical of the 
dimensionless bulk temperature parameter which has been displayed in Eq. (23). 
From this it follows that the Nusselt number depends only on the spacing ratio 
s/r 0 . Utilizing the numerical values obtained from evaluating Eq. (22), a 
graph of Nu^ as a function of s/r Q has been prepared and is presented as 
the solid curve of Fig. 3 - left hand ordinate. It is seen from the figure 
that Nu^ monotonically increases as s/r Q decreases^ and with this, it fol- 
lows from Eq. (25) that the wall-to-bulk temperature difference is smaller at 
smaller spacings for a fixed heat input. This trend is in accord with physical 
reasoning. 

Now, with the relation between heat transfer and temperature difference 
established by Fig. 3, we may return to the computation of the wall temperature 
variation as given by Eq. (20). Dividing through by Q/k and introducing 
Eq. (25), there is obtained 

*w - T° _ z 1 

Q7k “ (w/|i ) Pr + ir Nu^ 


( 26 ) 
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At a fixed flow, the value of (T w - T 0 )/(Q/k) at a particular z will diminish 
as the spacing diminishes. But, there will, of course, he a price in pressure 
drop to he paid to maintain the fixed flow. On the other hand, at a fixed pres- 
sure drop, w will diminish with decreasing spacing; and so, under this con- 
dition, the two terms of Eq. (26) vary in opposite directions. The dependence 
of the flow rate on spacing ratio is more acute (see Fig. 6, ref. l) than is 
that of the Nusselt number. So, for z values sufficiently large to establish 
fully-developed conditions, it would he expected that (T w - T 0 )/(Q/k) would 
increase with decreasing spacing when the pressure drop is fixed. 

Before leaving this section, it may he worthwhile to discuss briefly two 
other aspects of the results. First, it may he of interest to compare the heat 
transfer performance with that for laminar flow inside a circular tube 
(Nu^ = 48/ll). If we ask for the same wall-to-bulk temperature difference at 
a given heat flux per unit length, then this is achieved at a spacing ratio of 
2. 03. The second topic relates to the equivalent diameter. It is reasonably 
well-accepted that for laminar flow, the use of the equivalent diameter as a 
characteristic dimension is not sufficient to remove the dependence of the heat 
transfer results upon geometrical parameters. However, in the present problem, 
we are in a position to once again check on this point. From the definition 


4( cross-sectional area) 
heated perimeter 


(27a) 


it is easy to show that 


<*e _ 6( s / r 0 ) 2 
d V3 n 


(27b) 


Then, since Nu^ = (Nu^)(d e /d), we may find the dependence of Nu^ upon 
e 0 

s/r 0 by using the results which have previously been presented. A curve de- 


rived in this way has been plotted as a solid line in Fig. 3 - right hand 


ordinate. It is seen that the geometrical factors involved in the equivalent 
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diameter have given Nu^ a completely opposite trend to that of Nu^. How- 
ever, the variation in Nu^ is not less than that experienced hy Nu^. 

Peripheral variation of local heat transfer . - In the previous section, 
attention has teen focused on overall quantities such as the heat transfer per 

unit length, the hulk temperature and the Nusselt number. Here, consideration 
is given to the 

local heat transfer at various positions around the periphery of a cylinder. 

Since the velocity distribution varies with angular position around a cylinder 
while the surface temperature is prescribed to be independent of angle, it is 
evident that there may be a circumferential variation of the local heat transfer. 
To compute the local heat transfer, we use Fourier's law 



Utilizing the tabulated values of cjdj from table I with the and M from 

Ref. 1, the circumferential variation of q has been evaluated. The results 
are plotted on Fig. 4 as a function of 0 for parametric values of the spacing 
ratio s / r Q . By the symmetry considerations already mentioned in connection 
with Figs. 2(a) and (b), the variation over the entire circumference is known 
once the results for a typical 30° section have been given. 
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From Fig. 4, it is seen that for close spacings (small s/r Q ), there is a 
significant circumferential, heat transfer variation; while for open spacings 
(large s/r 0 ), the heat transfer is essentially uniform around the circumference. 
Clearly, for all spacing ratios of 2 or greater, there is essentially no dis- 
tinction between the case of peripherally uniform wall temperature or periph- 
erally uniform heat transfer. For many applications, this statement might he 
considered applicable even to spacing ratios as low as 1. 5. The trend of de- 
creasing peripheral variations with increasing spacing may be explained on 
physical grounds. For close spacings, a given cylinder is much influenced by 
the presence of its neighbors, with the result that the velocity distribution 
is circumferentially nonuniform and so is the heat transfer. As the spacing 
grows larger, the flow about a given cylinder is less disturbed by its neighbors. 
The resulting trend of the velocity toward axial symmetry is accompanied by a 
trend toward a more uniform surface heat transfer. 

It is interesting to note that the largest heat transfer occurs at 9 = 30°, 
the location of the largest open area for flow. On the other hand, the smallest 
heat transfer occurs at 9 = 0°, which is the location of the smallest open area 
for flow. From this, it may be inferred that if the heat transfer were pre- 
scribed to be circumferentially uniform while the surface temperature is per- 
mitted to vary, the hot spot would occur at 9 = 0°. 

SIMPLIFIED ANALYSIS FOR LARGE SPACINGS 

Both from physical reasoning and from the results of the foregoing analy- 
sis, it has been demonstrated that circumferential variations become smaller 
and smaller as the spacing ratio s/r Q increases. With this in mind, it seems 
worthwhile to reanalyze the problem using a model which neglects circumferential 
variations; and then, by comparing the results with those of the complete formu- 
lation already presented, to find for what conditions the simplified analysis 


is valid. 



E-725 


- 16 - 


A flow model having the property of angular symmetry with respect to a 
cylinder is an annulus surrounding the cylinder. The inner radius of the annu- 
lus is r Q , while the outer radius r* is chosen so that the flow area is 
identical to that in the actual configuration. From Fig. 2(a), it is seen that 
the actual flow area associated with a given cylinder is composed of 12 of the 
blackened typical elements. So, equating areas as follows 

*(r* 2 - r §) = 12^- tan £ - * jfJ 

there is obtained 

r * . (^.f .. 


(30a) 


(30b) 


The velocity problem for such an annulus is described by 

S 2 u + 1 Su _ 1 dp 
£ r 2 r Sr p, dz 

u(r Q ) =0, (r*) = 0 


(31a) 

(31b) 


It is easily verified by direct substitution that the solution for the velocity 
distribution is given by 



(32) 


which is identical to the nonseries part of the more complete velocity solu- 
tion (3). From this, the flow parameter M which has been defined in Eq- (4) 
can be evaluated to be 


M 1 i r 
M = — In — 
dit r r 


. jl£° V 

12 \ s / 96^ s/ 


8* T 19\ o / _ QRl a / ( 33 ) 

Next, by modifying the governing Eq. (5) to apply to the annulus, the 
temperature distribution problem may be formulated as 


Sr 2 + r 5? 


Ms* 


to -n - i 

A) 2 . (A) 2 ]] 

o 

u 

\sj Vs/ JJ 


(34a) 


T(r 0 ) = T v , || (r*) = 0 


(34b) 
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From this, the solution for the temperature is found to he 



It may he seen hy inspection that this result is identical to the nonseries 
portion of the more complete temperature solution (15). 

To compare the overall heat transfer characteristics with those previously 

obtained, we compute the wall- to -bull: temperature difference. Applying the 

definition of the bulk, temperature from Eq, (2l) to the annulus, there is 

obtained, /»r* 

2rtp J (T w - T)ur dr 

Tw - Tt = — - (36) 

By substituting from Eqs. (32) and (35), the integrand can he evaluated and 
then the integration carried out. The same type of functional relation which 
has been noted in Eq. (23) has also been found here. The end result of the 
calculation is reproduced in the Appendix to facilitate numerical calculations 
for s/r Q values which lie outside the range reported here. 

With the wall- 1 o -bulk temperature difference at our disposal, the Nusselt 
number Nu^ can be evaluated from Eq. (25). The results thus obtained have 
been plotted as a dotted line on Fig. 3 - left hand ordinate. By comparing 
with the solid curve from the more exact calculation, it is seen that the neg- 
lect of circumferential variations leads to a smaller wall-to-bulk. temperature 
difference at a fixed heat transfer per unit length. So, for a given bulk 
temperature at a cross section, the wall temperature prediction from the sim- 
plified calculation is low. If errors up to 5 percent in the Nusselt number 
prediction can be tolerated, then from Fig. 3 it is found that the results 
based on the neglect of circumferential variations may be used for spacing 
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ratios as low as 1.5- As the spacing becomes progressively smaller , the devia- 
tions become increasingly larger. The same comments apply to the Nu^ curves 
associated with the right hand ordinate. 

An analysis using the equivalent annulus approximation has been carried out 
in Ref. 3 for turbulent flow for spacings down to 1. 375. If the flow is fully 
turbulent, then the effects of circumferential variations will likely be smaller 
than in a laminar flow. On this basis, we may conclude that the use of the 
equivalent annulus model appears justifiable in the turbulent case for the 
spacings considered in Ref. 3. 


TABLE I. - LISTING OF o>j VALUES 


s / r o 

oijXlO 2 

oogXlO 3 

003 XIO 3 

< 04 X 10 ^ 

4.0 

2. 9649 

0. 39239 

0. 579 

~ 0 . 3 

2.0 

1. 2712 

. 12031 

.136 

-.05 

1.5 

. 71081 

. 014556 

-.195 

—1.5 

1.2 

. 37039 

-. 1323 

- 1 . 20 

— 6.0 

1.1 

. 25961 

-. 2065 

- 1. 12 

- 7.0 
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APPENDIX 

(a) Integral table. 
All integrals are over the range 0 to n/6 


/ 
f 

f 
/ 

f e)°°g MM = 0.01207813 

A 
A 
A 
A 
A 
A 


0 tan 0 ae = 0.05069781 

(H. cos 9)cos 68 le _ 0-0370858 
(cos 0) 8 

(in cos 6 ) COS 68 . dfl _ 0 . 04655214 
(cos 8) 10 

(in cos 8) cos 69 d8 _ 0 . 00966175 
(cos 0)- 4 

os 0)co 
(cos 0) 

cos 60 d0 = 0.195504853059 


= 0.0487139289629 


= -0.278775438621 


= 0.0270632938683 


= -0.3833738263351 


(cos 

0) 

8 

cos 

60 

d0 

(cos 

0) 

-4 

cos 

60 

d0 

(cos 

0) 

10 

cos 

60 

d0 

(cos 

0) 

-2 

cos 

60 

d0 

(cos 

0) 

12 " 

cos 

120 

d0 

(cos 

0) 

-6 

'cos 

120 

d0 


= -0.01087364485 


= 0.02483882327 


S (cos 0) 4 
J* (cos 0) 10 d0 = 0.354917462776 

J' (cos 0)“ 16 


/ 

I 


I 

f 

/ 

/ 

/ 

/ 

/ 


(In 

cos 

0) 

cos 

60 d0 : 


cos 

0) 

cos 

120 

d0 


(cos 

0) 14 



cos 

0) 

cos 

120 

d0 


( cos 

0) 16 


(In 

cos 

0) 

cos 

120 

d0 


(cos 

0)- 10 


12* 

cos 

0) 

cos 

120 

d0 


= -0.0434739 


= -0.0577564 


(cos 0) -8 
cos 120 d0 


= 0.000413828 


(cos 

0) ; 

L4 

cos 

120 

d0 

(cos 

0) 

-10 

cos 

120 

d0 

(cos 

0) 16 

cos 

120 

d0 

(cos 

0)' 

-8 

cos 

120 

d0 

(cos 

0) 18 

cos 

120 

d0 


(cos 0) 
8 


= 0.24953325973 


= 0.00934144946 


= 0.35013873273 


= -0.01065617196 


= 0.4834095306 


= 0.00993262699 


J ' (cos 0} 
/*( cos 0) 


d0 = 0.37912963362 
" 14 d0 = 1. 2290323215 


d0 = 1. 43544637803 
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Figure 1 


Flow configuration 
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(a) Cross-sectional view. 


(b) Typical element. 


Figure 2. - Diagram used in mathematical formulation. 
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Figure 3. - Nusselt number results. 
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Figure 4. - Peripheral variation of vail heat transfer. 
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